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We consider toy cosmological models in which a classical, homogeneous, spinor field provides 
a dominant or sub-dominant contribution to the energy-momentum tensor of a flat Friedmann- 
Robertson- Walker universe. We find that, if such a field were to exist, appropriate choices of the 
spinor self-interaction would generate a rich variety of behaviors, quite different from their widely 
studied scalar field counterparts. We first discuss solutions that incorporate a stage of cosmic 
inflation and estimate the primordial spectrum of density perturbations seeded during such a stage. 
Inflation driven by a spinor field turns out to be unappealing as it leads to a blue spectrum of 
perturbations and requires considerable fine-tuning of parameters. We next find that, for simple, 
quartic spinor self-interactions, non-singular cyclic cosmologies exist with reasonable parameter 
choices. These solutions might eventually be incorporated into a successful past- and future-eternal 
cosmological model free of singularities. In an Appendix, we discuss the classical treatment of 
spinors and argue that certain quantum systems might be approximated in terms of such fields. 
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I. INTRODUCTION 

At least since the advent of the first inflationary mod- 
els p|, cosmologies containing classical scalar fields have 
received widespread attention in the literature. From a 
purely phenomcnological point of view, such scalar fields 
are general enough to accommodate a rich variety of be- 
haviors. From a theoretical point of view, their invariable 
appearance in various theories of nature makes them nat- 
ural candidates for cosmological applications. In spite of 
these facts, one might wonder to what extent scalar fields 
are singled out by the former considerations. Could other 
classical homogeneous fields play a significant role in cos- 
mology? 

In the present paper, we consider the possibility that 
classical, homogeneous, spinor fields might play a role 
in cosmology. By a classical spinor field, we simply 
mean a set of four complex-valued spacetime functions 
that transform according to the spinor representation 
of the Lorentz group. Although the existence of spin- 
1/2 fermions is both theoretically and experimentally 
undisputed, these are described by quantum spinor fields. 
It is unclear when fermionic quantum fields might be 
consistently treated as classical spinors. It is generally 
held that there exists no classical limit for fundamental 
quantum Fermi fields; however, one can imagine clas- 
sical spinors as arising from an effective description of 
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a more complex quantum system. We address possible 
justifications for the existence of classical spinors in an 
Appendix. For the bulk of this paper, we will simply 
presuppose their existence. 

We find that classical spinors are mathematically per- 
fectly consistent. Physically, one might object that 
spinors violate Lorentz invariance and isotropy. With- 
out being explicit about the precise nature of Lorentz 
transformations in a general gravitational background, 
let us point out that Lorentz invariance is broken in 
any Friedmann-Robertson- Walker cosmology, regardless 
of whether a spinor has a non- vanishing expectation value 
or not. On the other hand, we shall see that eventual vi- 
olations of isotropy caused by the spinor do not prevent 
consistent solutions of Einstein's equations, and might 
actually remain undetectable. 

Compared to scalar fields, spinor fields have attracted 
little attention in cosmology. One of the first papers 
about the subject was Taub's study of the Dirac equa- 
tion in various cosmological spaces 2] . Brill and Wheeler 
dealt with neutrinos in gravitational fields 3] . The quan- 
tization of a spinor field in an expanding universe was 
considered by Parker Q, and the quantization of grav- 
ity coupled to a spinor was addressed in [5|, |6j . Explicit 
solutions of the Dirac equation in an open Friedmann- 
Robertson- Walker spacetime have been considered in • 
Solutions of the Einstein equations coupled to a spinor 
in Bianchi Type I spaces have been extensively studied 
by Saha and Shikin 8]. 

The structure of this paper is the following. In Sec- 
tion^ we discuss how to couple a spinor to gravity. The 
reader familiar with the formalism might want to skip 
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to the next section and eventually refer back for nota- 
tional details. Section 11111 deals with the basic cosmo- 
logical equations and general solutions in terms of an 
arbitrary spinor self-interaction term. In Section IIVI we 
study inflation driven by a spinor field and compute the 
spectrum of primordial density perturbations. Section 
Section [V] presents a cyclic, non-singular model of the 
universe that critically relies on the properties of a spinor 
field. Finally, in Section lVll we summarize our results and 
draw our conclusions. The Appendix comments on the 
meaning and properties of classical spinors. 



II. FORMALISM 

In this section we briefly review how a spinor field 
is coupled to gravity. For complete discussions about 
spinors in curved spacetimes, see [9l lTollTl| . 

Because the group of diffeomorphisms does not ad- 
mit spinor representations, in order to couple a spinor 
to gravitation one introduces the Lorentz group (which 
does actually have spinor representations) as a local 
symmetry group of the theory. Under diffeomorphisms 
x n x^(x u ), a spinor ip is a scalar, i/; — > t/> = ip, but 
under a local Lorentz transformation with parameters 
Xab(x) a spinor transforms according to 



that preserve the orthonormality relation © at each 
point. Thus, the spacetime metric only determines the 
vierbein up to such local Lorentz transformations. For 
this reason, one must ensure that any Lagrangian formed 
with the aid of the vierbein is invariant under the Lorentz 
group acting as a local symmetry. Invariant terms con- 
taining derivatives of a spinor can be constructed through 
the covariant derivative 



(5) 



which transforms as a (covariant) vector under diffeomor- 
phisms and as a spinor under local Lorentz transforma- 
tions. The 4x4 matrix f2^ is the spin connection 



a. 



i 



-W^afeE , UJfiab = aS ^vb ■ 



(6) 



The coefficients LJ^ab arc the Ricci rotation (or spin) co- 
efficients. 

The vierbein and the flat-space gamma matrices allow 
one to define a new set of gamma matrices 



(7) 



that satisfy 1 the algebra {Ti i ,T v } = 2g» v . These can be 
used to write down a generalization of the Dirac action 
in a curved spacetime background, 



ip — > ip = exp 



-A a6 (a;)S 



a b 



(1) 



d x e 



D^r^) - v 



(8) 



where S ab = |;[7 a ,7 fc ] are the generators of the spinor which we have written in a symmetrized form. Here, e 
representation of the Lorentz group, and the 4x4 ma- 
trices 7 satisfy the Clifford algebra {j a ,"f b } = 2r/ ab . We 
shall choose the Dirac-Pauli representation 



7° = 



1 
-1 



7 = 



a 1 
-a 1 



(2) 



where the a 1 are the conventional 2x2 Pauli matri- 
ces. Then, 7 = (7°)^ is Hermitean, and the 7* = 
— (7*)' are anti-Hermitean. For later convenience, we 
shall define the additional (Hermitean) gamma matrix 

75 = «7o7i7273- 

A fermion is coupled to gravitation with the aid of a 
"vierbein" e M a , a set of four contravariant vector fields 
that satisfy the orthonormality condition 



9fj,v^- a& b — 'Hab ■ 



(3) 



where g M „ is the spacetime metric and rj a b is the 
Minkowski metric f] a b = diag(l, — 1, — 1, — 1). Latin in- 
dices enumerate each of the vectors in the vierbein while 
Greek indices enumerate the spacetime components of 
each of these vectors. Spacetime and Lorentz indices are 
raised and lowered with the spacetime and Minkowski 
metrics respectively, leading to associated sets of vectors 
such as e afl and e M °. 

Local Lorentz transformations A(x) are just local 
"reshufflings" of the vierbein vectors 



is the determinant of the vierbein e M a , and the Dirac ad- 
joint %[) is given by ip^° . The covariant derivative acting 
on the adjoint is D^tp = d^ij) — V'^m- By an integration 
by parts, the kinetic term of the spinor can be cast in the 
"conventional" form iipT^D^ip. The term V stands for 
any scalar function of ip, tp and possibly additional mat- 
ter fields. When a particular form of V is later needed, 
we will assume that V only depends on the scalar bilinear 
ipip. It turns out that this choice is general enough for 
our purposes. More general interactions in Bianchi Type 
I spacetimes have been considered in the series of papers 


The Lagrangian iJSJ describes how the spinor is cou- 
pled to gravity, but it does not specify the dynamics of 
gravity. We shall assume that the latter is governed by 
the Einstein-Hilbert action. Hence, we consider a spinor 
minimally coupled to general relativity, 



S 



(9) 



where R is the scalar curvature, S^p is given by Eq. 
and S m describes additional matter fields, such as scalar 



A, 



b p fJ, h 



(4) 



1 Note that, while the 7 a 's do not transform under local Lorentz 
transformations or diffeomorphisms, the T^'s do. They inherit 
their transformation properties from the vierbein. 
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fields or gauge fields. The symmetries we have postulated 
up to now, diffcomorphism and local Lorentz invariance, 
certainly allow for the presence of additional terms in 
the action. For example, we could have written down a 
non-minimal term like iptpR. However, as we are going 
to see, in an expanding universe ipip decays at least as 
fast as 1/a 3 . Therefore, during cosmic expansion such a 
term would quickly become negligible. This situation is 
in sharp contrast with the case of a scalar field <fi, where 
a priori there is no reason to expect a term proportional 
to 4>R to be negligible (see however 0,0). Eventually 
this fact could be relevant in models where a spinor field 
drives late time cosmic acceleration. 

Varying the action © with respect to the vierbein e p a 
leads to the Einstein equations 

G^u = 3T^ V , (10) 

where the energy momentum tensor is given by the 
variation of the matter action, 



For these isotropic solutions of the Einstein equations 
to exist, the energy-momentum tensor of the spinor must 
be compatible with the symmetries of the metric <|15[1 . 
homogeneity and isotropy. At this point, note that ho- 
mogeneity of a spinor is not a gauge-invariant concept; 
by a local Lorentz transformation (Q, it is always pos- 
sible to transform a homogeneous (space-independent) 
spinor ip(t) into an inhomogeneous (space-dependent) 
one ip(t,x). We are going to look for spinor solutions 
of the Dirac equation that can be written as a gauge- 
transformed homogeneous spinor. If that is the case, 
there exists a vierbein where the Dirac equation allows 
space-independent solutions. In our case such a vierbein 
is given by 2 



e" = « M o ,ef 



1 



(16) 



In the gauge l|16|l the equation of motion of a space- 
independent spinor (|14J) reads 



T — 



e 5e v a 



(11) 



Note that we work in units where 8irG/3 = h = c = 1. 
Substituting the action into Eq. I jlljl we obtain, after 
an integration by parts, the energy momentum tensor of 
the spinor (on-shell), 

W Tw = \ [#VA,)V> - D^T^] - g^L 4 , . (12) 

We have used a relation that follows from the Lorentz 
invariance of the spinor Lagrangian, 

(13) 

to rewrite the result of the spinor variation. On the other 
hand, varying the action with respect to the held ip yields 
the equation of motion of the spinor, a generalization of 
the Dirac equation to a curved spacctimc, 



dV 



0. 



(14) 



If the action is real, the variation of the action with re- 
spect to ip yields the adjoint of the previous equation. 



III. COSMOLOGICAL SOLUTIONS 

Because we are interested in cosmology, in this pa- 
per we deal with homogeneous and isotropic FRW space- 
times. Current observations favor a flat universe [lj], so 
we assume the spacetime metric to be spatially flat, 



ip+-Hrp + vr°V'il) = 0, 



(17) 



where a dot (') denotes a time derivative, a prime (') 
denotes a derivative with respect to iptjj, and H = 
d(\oga)/dt is the Hubble parameter. The equation man- 
ifestly admits space-independent solutions, and hence, 
spinor observables like the energy momentum tensor are 
also homogeneous. 

One should also verify whether spinors are compatible 
with the isotropy of the FRW metric. The q 1 Einstein 
Eqs. = Gq 1 — To 1 are satisfied only if Tq % vanishes. 
This is possible for conventional matter forms (perfect 
fluids and homogeneous scalars), but it is not generally 
true for a spinor. In fact, in spatially open or closed uni- 
verses, it is not possible to satisfy the constraint Tq 1 = 
unless 4>ip is zero 0- Q In a spatially flat universe how- 
ever, the equation of motion (|17f) automatically implies 
the vanishing of To 1 , so that the presence of the spinor is 
consistent with the isotropy of the metric. 

A convenient combination of the remaining Einstein 
Eqs. i|l(J|) . the o° and the ^ , involves the energy density 
Pk and the pressure pk of the different constituents of the 
universe, 



H 2 = p^ 



'-.\ftiP + Pm+ i{P4, + Pm)}a. 



(18) 
(19) 



Here, the sub-index ip stands for the spinor and m for 
any additional matter, such as dust, radiation, or even 
dark energy. The spinor's energy density and pressure 



ds 2 



dt 2 - a 2 (t)dx 2 



(15) 



2 In open or closed FRW universes, one can construct a vierbein 
that allows homogeneous spinor solutions. These are formed 
from the Killing vectors of the homogeneous 3-dimensional spaces 
of Bianchi Type V and Type IX, respectively. 
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are given by the corresponding components of the energy 
momentum tensor (I12I1 . 



V. 



(20) 
(21) 



The equation of state of the spinor is the ratio of its 
pressure to energy density, and hence, it is given by 



w 4: = — = 77 ■ 

Pi, V 



(22) 



The equation of state is not restricted to be within the in- 
terval — 1 < w < 1. For a conventional massive fermion, 
V = rmpijj, the equation of state agrees with that of a 
fluid of dust, = 0. For more general choices of V, 
may acquire any real value. 

It is possible to directly integrate the spinor equation 
of motion if V only depends on ipip. It follows directly 
from the Dirac equation i|17[l that 



A 

73- 



(23) 



where A is a time-independent constant. Note that this 
result is valid for any time dependence of the background 
geometry, a(t), and thus, is valid regardless of the dom- 
inant energy component of the universe. In an expand- 
ing universe, the value of ipip monotonically decreases, 
whereas in a contracting universe the value of iptj) mono- 
tonically increases. These facts do not imply however 
that the energy density of the spinor follows the same 
behavior. The Dirac equation can be cast as a continuity 
equation 



p + 3Hp(l 



0. 



(24) 



Integrating Eq. I|24|) or directly from Eq. I|23l) it is pos- 
sible to find p4, as a function of the scale factor a for 
arbitrarily given Vtyip), 



Pi> 



V 



ifjip—A/a 3 



(25) 



Conversely, given an arbitrary function p^ia) one can 
always find a Vfyip) such that p(a) is a solution of the 
equation of motion l|24(l . 



V(ipip) = Pi, 



a=(A/i>i>) 1 /f 



(26) 



In conclusion, a spinor field can accommodate any de- 
sired behavior of its energy density by an appropri- 
ate choice of V. In that respect, a spinor field is 
completely different from a scalar field. A (canonical) 
scalar field cannot violate the null 3 energy condition ^(| , 
whereas a spinor field can violate any desired — weak 4 , 



null, strong 5 or dominant 6 — one. Even non-canonical 
scalar fields — "k-fields" 01 — cannot reproduce the be- 
havior of a spinor. In the former there are barriers that 
prevent a transition from / o + p>0top + p<0, whereas 
such barriers are nonexistent for spinors. The converse is 
however true. A spinor can reproduce the behavior of a 
scalar field. In particular, it can drive inflation and late 
time cosmic acceleration. 

Although these solutions of Einstein equations sourced 
by a spinor are perfectly valid and consistent, they might 
break isotropy. By that we mean that the spatial compo- 
nents of certain vector quantities that involve the spinor 
do not necessarily vanish, and hence, are not invariant 
under spatial rotations. For instance, it turns out that 
for non-trivial homogeneous solutions of the Dirac Eq. 
(|17|l the spatial components of the vector j* 1 = ijjT^t/j are 
generically non-zero. If the action © does not include 
a coupling between j M and any other observable vector 
quantity, such a violation would be undetectable. On 
the other hand, if the action contained such a coupling, 
there still exist some spinors for which j % = 0, such as, 
for instance, 



^=(^1, 0,0,0), 



(27) 



Note that this form of the spinor is compatible with the 
equation of motion Q17[l. In other cases, no choice of 
a spinor prevents isotropy violations. There is no non- 
trivial spinor such that the "pseudo-vector" {jj^^^ij} has 
non- vanishing spatial components 0, @. But again, if 
there is no coupling between the latter vector and any 
other observable component (say, because of parity con- 
servation), such a violation would remain undetectable. 

Although in this section we have mainly assumed that 
V only depends on the scalar bilinear ijnj), some of the re- 
sults can be easily generalized for rather arbitrary choices 
of V. Consider, for example, any V that is invariant un- 
der the global transformation if) — > e la ^j 1 for arbitrary 
constant a. Such a symmetry means that the fla- 
vor is conserved, and hence, there is a conserved current 
V ^(ipT^ip) = 0. Then, for a homogeneous spinor 



A 

"07V = -r, 



(28) 



which already suggests that Eq. I|23() is not just a con- 
sequence of our choice of V. In fact, writing down the 
4-spinor in terms of two 2-spinors, ip — (u,v), it follows 
from the identity 

(^V) 2 + (#75V>) 2 + (^7°75^) 2 = 
(Vn-V) 2 - 4 [(v!u)(v*v) - (u^v)(v^u)] 

and the Cauchy-Schwarz inequality that 



(29) 



(V^) 2 + (#75V-) 2 + (^7°75^) 2 < ^r. (30) 



V > o 

- p > and p > 



' p + p > and p + 3p > 
' P>\P\ 
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Therefore, spinor bilinears generically decay during the 
expansion of the universe, without regard to the precise 
form of the spinor interaction. In particular, in an ex- 
panding universe, there are no non-trivial solutions of 
the spinor equations of motion with constant tp. 



IV. INFLATION 

In this section, we investigate the possibility that a 
classical spinor field could drive inflation. A sufficiently 
long stage of inflation p| explains many of the features 
of our universe that remain unexplained otherwise (see 
[Tfij for diverging claims). Nevertheless, the nature of 
the component that was responsible for inflation remains 
unknown. Most inflationary scenarios rely on a homoge- 
neous scalar field rolling down an appropriate potential; 
however, at present there is no direct experimental ev- 
idence for the existence of fundamental scalar fields in 
nature. Hence, a natural question is whether a another 
type of field could have driven a stage of inflation in the 
early universe. Inflation driven by a vector field has been 
considered by Ford |23|. and inflation driven by a spin- 
ning fluid has been discussed by Obukhov pof. 



A. Background 

By definition, inflation is a stage of accelerated expan- 
sion of the universe, a > 0. It follows from Eq. i|19|) 
that any component driving inflation has an equation of 
state that obeys w < —1/3. (We assume p to be posi- 
tive.) Three types of inflation are mainly considered in 
the literature: pole- like inflation with w < —1, de Sit- 
ter inflation with w = — 1 and power- law inflation with 
— 1 < w < —1/3. It is easy to verify from Eq. 11' I'll that 
inflation (or simply expansion) with constant equation of 
state w results from a "potential" 



V = (W) 



l+w 



(31) 



If w < —1 the expansion runs into a future singularity, 
while if w > — 1 the expansion runs into a past singular- 
ity. However, by an appropriate choice of V, the universe 
could pole-like inflate in the past and power-like inflate 
in the future. 

If w — —1, the formula l|31|l implies that V is constant, 
as for a cosmological term. It is possible to relax the con- 
dition on the function V by looking for a stage of nearly 
de Sitter inflation, w « — 1. In terms of the function V, 
the condition for nearly de Sitter inflation is 



dlogV 



d log ipip 



< 1. 



(32) 



Note that the latter condition alone suffices to guarantee 
quasi de Sitter inflation. This is to be compared with con- 
ventional scalar-field driven inflation, where two slow-roll 
conditions are needed. In general, any V that asymptotes 



to a positive constant at large 'tjj ip satisfies Eq. (|32ll . Ex- 
amples of such functions V are log[l + (-0VO™L tavii n-rpip 
and ("ipip) 11 /(I + V'VO™ f° r arbitrary positive n and suffi- 
ciently large iptp. 

Although nothing prevents a spinor field from driving 
inflation, certain facts make this possibility unappealing. 
Inflation solves the homogeneity problem if it lasts for 
about 60 e-foldings. Let us assume that V is such that 
w < —1/3 for tptp > (if)ip) e and w = —1/3 for iptf) = 
('tpip) e (Fig. The end of inflation is determined by 
(ipip) e , the point where the equation of state w crosses the 
^critical" value -1/3. For instance, for V = {$ij]) n /{l + 
tpip) n inflation ends once ip-ip reaches 3n/2 — 1. If the 
initial value of the scalar bilinear is (iptp)i, then inflation 
lasts a number of e-foldings N given by 



AT 1 1 

N= * log m:- 



(33) 



It follows that during 60 e-foldings -t/jip changes by eighty 
orders of magnitude! This is to be compared with a con- 
ventional chaotic model, where the scalar field changes by 
just an order of magnitude. This fact is particularly im- 
portant in nearly de Sitter inflation, since the "flatness" 
condition l|32|l has to be satisfied for a range of values of 
t/jip that encompasses eighty orders of magnitude. 

An important difference between inflation driven by 
a spinor and the conventional scenarios is the reheating 
mechanism after the end of inflation. In the conventional 
scenarios, the universe is reheated when the scalar field 
starts oscillating around the bottom of its potential and 
decays into particles 0, |^, • If inflation is driven by 
a spinor field, the quantity tptp evolves according to Eq. 
(|23|l and hence does not oscillate. Nevertheless, there are 
several mechanisms to reheat the universe. One of them 
is gravitational particle production at the end of inflation 
[23l |24| ; more efficient ways have been suggested in [25| . 



B. Perturbations 

One of the most appealing features of many of the 
conventional inflationary models is their prediction of an 
adiabatic, nearly scale invariant spectrum of primordial 
density perturbations, in agreement with current obser- 
vations. Our goal in this section is to compute the power 
spectrum of density perturbations generated during a 
stage of nearly de Sitter inflation driven by the spinor 
field f/;. In a proper treatment of the problem, we would 
perturb both metric and spinor and solve the linearized 
Einstein equations. The nature of the spinor makes this 
path cumbersome, so we shall rely on a simplified anal- 
ysis, where we only perturb the spinor in a given, fixed, 
spacetime background (de Sitter spacetime). 

We shall characterize density perturbations Sp by the 
variable 



Sp 



(34) 
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FIG. 1: A plot of a generic interaction that yields inflation. 
The corresponding equation of state is also shown in the dia- 
gram. For large values of ipip, the interaction is flat, Eq. 132L 
allowing nearly de Sitter inflation . At the critical value (ipip) 
the equation of states reaches — 1/3 and inflation ceases to be 
possible. 



This quantity is somewhat analogous to the Bardeen vari- 
able, which is conserved on large scales in the absence of 
entropy perturbations, and which can be directly related 
to the cosmic microwave background temperature fluc- 
tuations. The source of the density perturbations Sp are 
the fluctuations Sip of the spinor field around its homo- 
geneous background value ipo. We treat ipo as a classical 
field, and the fluctuations Sip as a quantum field in an 
expanding universe 

H = ^372 j d 3 k ]T [u(t, k, a)a(k, oY ls + 



— iki 



(35) 



The index a runs over the two spin states of the 
spinor, and the operators a and b are particle and an- 
tiparticle annihilation operators, {a(k, a), a^(k', er')} = 
{b(k,a),tf(V,a')} = S^(k~ k')8 aa .. 

The power spectrum V(k) is a measure of the fluctu- 
ations of the variable £ on comoving scales of size 1/fc, 
and it is implicitly defined by the relation |26j 



(((t,x)((t,x + r)) = 



dk sin kr 
k kr 



P(k). (36) 



Here, ( ) denotes the expectation value in an appropri- 
ately chosen vacuum state, a|0) = fe|0) = 0. Using ex- 
pressions (I20|) and l|21|) for the energy density and pres- 
sure of the spinor field respectively, we find that 



c 



Sip ip + ip Sip 
ipip 



(37) 



where we have dropped the subindex that denotes back- 
ground quantities. Substituting Eq. lj3"Tj) into the left 
hand side of Eq. Ip^fljl we obtain (for r = 0) 



(C(t,x)((t,x)) = 



2(Sip(t,x)ip(t) ■ iP(t)Sip(t,x)) 



(38) 



where we have used the fact than only terms with equal 
number of creation/annihilation operators have a non- 
vanishing expectation value. 

Using the Pauli-Fierz identities [27j we can express the 
previous four spinor expectation value in terms of pertur- 
bation bilincars, 



(CC> = 



(SipSip) (ip'ya*P)(Sip j a Sip) 



2ipip 



2(#) J 



(39) 



Note that the second term in the right hand side intro- 
duces violations of isotropy in the power spectrum unless 
ip^f a ip vanishes 7 . Because we are only interested in an 
estimate of the amplitude and the k dependence of the 
spectrum, we can concentrate on the first term on the 
right hand side. Substituting the expansion 135f) into 
that term we finally obtain that the power spectrum is 
of the order 



V(k) 



k 3 y-v v(t, fc, a) v(t, fc, a) 



At: 2 



(1P1P) 



(40) 



The time evolution of v is dictated by the equation of 
motion of Sip. The field Sip itself satisfies the linearized 
Dirac equation iT°DoSi/j + iT l DiS^ — mSip = 0, where 
we assume that m = V' is small but non-zero and V" is 
negligible. It is convenient to work with the rescaled field, 
Sip = a 3 / 2 Sip instead of Sip. The rescaled field behaves as 
a spinor with a time-dependent mass in flat space, and 
in particular, v satisfies 



«7 



, dv 
dr) 



+ 'fkiv — ami' = 0, 



(41) 



where rj denotes conformal time, dr] 
e~ Ht /H runs from 



dt/a. In de Sit- 
ter space, rj = —e~ nt /H runs from — oo to 0, and 
a = —1/(Ht]). Solutions of the Dirac equation l|41|) in 
a de Sitter background were studied by Taub 2]. The 
ansatz v = (w+Vf, V-V-), where V+ and V- are two 
time-independent two-component spinors, yields the sec- 
ond order differential equation 



+ [fc 2 + a 2 m 2 ± i(am)']v± = 0. 



(42) 



Different linear combinations of the solutions to Eq. 142|) 
correspond to different choices of vacuum state. We 



7 The power spectrum is isotropic if the Fourier transform of the 
correlation function on the lhs of Eq. 1361 only depends on 
k = [fc| , and not on k itself. For simplicity, we have implic- 
itly assumed isotropy in our definition of the power spectrum V, 
m the rhs of Eq. 
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choose the standard prescription where v± cx e lkr> as 
T] — > — oo The corresponding properly normalized 

spinor solutions are then 



-nkr) e ^ m / 2H 



-irkr) e 



i/2H 



_ff< 2) (-fcr;)fe3/fc 
e *m/H H (2) { _ kri) 

Hl 2 \-kr 1 ){k 1 -ik 2 )/k" 
-Hl 2 ~>(-kT,)k 3 /k 

6 ™' H Hf(-t,) 



(43) 



The functions tf^ and #£ 2) arc the Hankcl functions of 
the first kind and second kind [2^], and v = | — im/H. 

Up to the factor a -3 / 2 , the previous spinors oscillate as 
e ikri £ or moc j es ingide the horizon, — krj ^> 1. Using the 
asymptotic expansion for the Hankcl function in the limit 
— krj <C 1 (modes outside the horizon) 



-TM 



-krj 



(44) 



it is straightforward to verify that the power spectrum 
"freezes" on large scales and becomes equal to 



V(k) 



k ITV M 2 • U 77,71 

2^4 smh ^" 



(for — krj <C 1). 

(45) 

Such a power spectrum has spectral index n = 4, in 
strong disagreement with experimental results consistent 
with a scale invariant spectrum with n « 1 14] . The 
constant A is the quantity that appears in Eq. (|23|1 . 
Equation l|45|) can also be used to estimate the power 
spectrum of the density contrast, 



(46) 



Because during inflation V is nearly constant while a 
grows exponentially, spinor fluctuations are highly sup- 
pressed with respect to, say, scalar field density fluctua- 
tions. 

The k 3 dependence of the power spectrum il4*5|) and 
the a 3 decay of the density contrast are to some extent 
an expression of the conformal triviality of the system. 
Indeed, a massless spinor is conformally invariant, and 
the power spectrum of a massless spinor in flat space- 
time displays the same fc 3 dependence. Our calculation 
shows that even the inclusion of a conformal symmetry- 
violating mass does not significantly alter this result. 
Note that although expression (|45|l vanishes for to = 0, 
this merely reflects the chiral asymmetry of the first term 
in the expansion (|39|l . In the limit of zero mass, the dis- 
carded terms give the dominant contributions, which also 
are proportional to k 3 . 

In conclusion, at the level of our simplified preliminary 
analysis, it seems that a stage of (quasi) de Sitter inflation 
driven by a spinor cannot seed a scale invariant spectrum 



of primordial density perturbations by itself. Eventually, 
a light scalar field present during inflation (as in curvaton 
models [l^) ma y solve this problem. 



V. NON-SINGULAR CYCLIC COSMOLOGIES 

One of the most intriguing issues in cosmology is the 
ultimate origin of the universe and the character of its 
initial state. One of the attractions of cyclic cosmolo- 
gies [2!| lMJ, |3l| is that — to the extent that they are truly 
cyclic, returning to the same state after each cycle — they 
dispense altogether with that problem. Since they are 
past eternal, there is no need to formulate initial con- 
ditions from which the universe is evolved into the fu- 
ture. Furthermore, the universe has always existed for 
the same reason, so that there is no need to ask where 
it originated from. However, many of the cyclic universe 
models that have been proposed so far suffer from sin- 
gularities that prevent a continuous account of cosmic 
history. At a certain time, the universe evolves into a 
singular state where the conventional low-energy effective 
theory description of the universe breaks down. Further- 
more, even if the singularity is regulated in some way, 
these models can still lead to inconsistencies [32|. 

In this section we describe a scenario which avoids this 
latter breakdown (see ^jj^ for alternatives). For simple 
choices of the self- interaction term VfyijS), cyclic cos- 
mologies free of singularities exist. Here, we present a 
simple model illustrating this point. Though we make no 
claim that this simple model leads to an entirely satisfac- 
tory cosmology, we see no obstacle to refining the basic 
idea into a more realistic description of our universe. 

Consider a spatially flat universe that contains "mat- 
ter" (dark energy, dust and radiation) and a homoge- 
neous spinor field ip, and suppose that the interaction 
term V in Eq. (JSJ has the form given in Figure [21 The 
"potential" is negative for "small" values of ipx/j, and it 
becomes negative and decreases "fast enough" for large 
values of ipip. Such an interaction might be given for 
instance by 



V(iptp) = A^p + mipip - A(V'V') 2 



(47) 



Here Ay, is a (negative) contribution to the total cosmo- 
logical constant, to is a (positive) mass and A is a (posi- 
tive) coupling constant. Hence, such a model describes a 
conventional, self-interacting, massive spinor with a neg- 
ative contribution to the vacuum energy. 

In order to describe cosmic evolution in such a uni- 
verse, let us arbitrarily start our description during ex- 
pansion. Suppose that the matter energy density domi- 
nates over the energy density of the spinor at a time when 
the latter is positive (region II in Fig. [21) In an expanding 
universe ipip ^ a ~ 3 is driven to values where the energy 
density of the spinor becomes negative, while the energy 
densities of radiation (cx a -4 ) and dust (oc a~ 3 ) are "red- 
shifted away" and tend to zero. The only assumption we 
have to make at this point is that the energy density of 
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V 




FIG. 2: Generic form of the potential V in a cyclic non- 
singular universe. 



dark energy does not increase at late times 8 . Then, if 
in Eq. I|47[) is large enough, there necessarily exists 
a value of the scale factor a max where the total energy 
density ptot becomes zero, ptot — (region I in Fig. [31) 
It follows from Eq. (|18(l that at a max , a — 0. In addition, 
at a max the right hand side of of Eq. 1)19(1 , 

a = 2 — Pm{w m -w^,), (48) 

is negative, since if ptot reaches zero from a positive value, 
the combination w m — has to be positive. Thus, at 
amax, a = and a < 0, so the universe automatically 
starts contracting. 

After the universe starts contracting, tpip oc a~ 3 re- 
verses its motion and starts growing. At the same time, 
the energy densities in matter increase as the universe 
contracts. At sufficiently small a we can assume that 
the matter component of the universe mainly consists of 
radiation. Assume that for large values of xptp, V de- 
creases faster than the rate at which the energy density 
of radiation increases, > 1/3 (recall that the universe 
contracts). If V is well approximated by a power at large 
values of -ipip, this implies, from Eq. I|31[l . that \V\ grows 
faster than (i^i/i) 4 / 3 , which is satisfied by the interaction 
(E7)l . Then, the ratio of spinor to radiation energy densi- 
ties steadily approaches —1 (region III in Fig. |3J) Again, 
there exists then a scale factor a m i n where ptot — 0. At 
&mim a = and from Eq. (|19fl a > since, when p tot = 



By a suitable modification of the interaction in Fig. [5] the spinor 
field could also account for dark energy. In that case one can 
drop the latter assumption and, up to the constraints on the 
parameters of our toy model, the rest of our discussion remains 
unaltered. 



and w > 1/3 are satisfied, Eq. i|19|) reads 

a « a min ^(3w^ - 1) > 0. (49) 

Hence, at a m i n the universe bounces and starts expand- 
ing until it again reaches region II. From that point on 
cosmic history as described above repeats itself. Note 
that cosmic evolution is singularity- free throughout. Af- 
ter the bounce from contraction to expansion, the total 
equation of state of the universe evolves from iu tot = — oo 
towards w to t = 1/3 during radiation domination. Hence, 
there is a stage of inflation between the bounce and radi- 
ation domination, and eventually a spectrum of nearly- 
scale invariant perturbations in the matter fields could 
be seeded. 

It is possible to estimate the values of the parame- 
ters A^,, m and A by imposing certain observational con- 
straints on our model. Let us set today's value of the 
scale factor to 1 and denote the present value of the spinor 
bilinear by (ipip)o- Requiring the energy density of the 
spinor to be subdominant today we find 

|A^ + m(M0o| « Pent « HT 121 , (50) 

where we have assumed that the term proportional to 
A is negligible today. The previous relation shows that 
our model requires a certain degree of fine tuning, but 
let us point us that the required fine tuning is of the 
same order as the one needed to explain late time cos- 
mic acceleration. If condition l|5U|) is satisfied, the spinor 
field remains subdominant all the way into the past, un- 
til the moment when its (negative) energy density ex- 
actly compensates the (positive) matter energy density 
and the universe bounces (Fig. |3J). Primordial nucle- 
osynthesis is the earliest epoch when conditions in the 
universe can be probed. In order for the bounce to occur 
before nucleosynthesis, and in order not to conflict with 
its standard predictions, the energy density of the spinor 
should be subdominant at that time as well. Since pri- 
mordial nucleosynthesis occurs at a ~ 10 -10 , and today 
(p r )o ~ W~ 4 p cr it, it follows that 

A(V^)o « 10~ 145 - (51) 

Finally, in order for our field-theoretic classical descrip- 
tion to remain valid throughout cosmic history, we im- 
pose that the energy density of radiation at the bounce 
be significantly below the Planckian energy density, 

A(^V)o » 10~ 188 - (52) 

Due to the freedom in the dynamical variable (ipip)o, 
there is a large set of parameters that satisfy the con- 
straints l(5Tj|) , lfST|l and ((52|) . It is easy to verify that a set 
of parameters that satisfies all the constraints is 

A w -1CT 2 • (l(T 3 eV) 4 , m » 10 _3 eV, A w l(T 5 GeV" 2 , 

(53) 

where we have assumed (i/'VOo = 10 95 , which implies 
that at the bounce (ipt/}) m in ~ 10~ 43 . Hence, throughout 
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log |p| 




FIG. 3: Schematic plot of the energy densities of matter (p m ) and spinor (p^) as a function of the scale factor a. In Region 
I the energy density of the spinor is negative. If the universe expands, it starts contracting at a m ax, when spinor and matter 
energy densities add to zero. In Region II the spinor remains subdominant. Once its energy density becomes negative again 
in Region III, its importance starts growing until the moment when the total energy density is zero again and the universe 
bounces back. 



cosmic history ijjtp remains much smaller than 1. Notice 
that the cosmological term in i|53|l is of the same order 
as the component that is presently driving cosmic accel- 
eration, the mass m agrees with common neutrino mass 
models, and the coupling constant A is of the order of 
Fermi's constant. From this point of view, the parame- 
ter choices l|53|l do not appear to be unphysical. 

Certainly, some aspects of our cyclic scenario still re- 
main to be discussed. We have not addressed the issue of 
how to obliterate the inhomogeneous debris of previous 
cycles, and we have not taken into account particle pro- 
duction at the bounce. Eventually, both issues might be 
related. Finally, let us point out that even in the absence 
of matter, the spinor can support a cyclic universe, where 
o-min and a max are simply determined by the values of iptp 
where V becomes zero. 



VI. SUMMARY AND CONCLUSIONS 

In a flat FRW universe there are consistent solutions 
of Einstein's equations coupled to a homogeneous spinor 
field. For such solutions, the scalar bilinear ^ip is propor- 
tional to a -3 , as for the number density of a gas of non- 
relativistic particles. The energy density of the spinor is 
given by an (a priori) arbitrary self-interaction term V. 
For a given form of the spinor energy density p^,(a), one 
can always find a self-interaction V{ipip) that has p-^(a) 
as a solution of the equations of motion. Thus, canonical, 
classical, homogeneous spinors can violate any desired en- 
ergy condition, and their behavior in general cannot be 
reproduced by a minimally coupled, homogeneous scalar 
field. 



A spinor field can also support a sufficiently long stage 
of inflation, provided the self-interaction term V satisfies 
a single condition on its slope for an exponentially large 
range of iptp. This condition is satisfied, for instance, 
if V asymptotes to a constant value at large values of 
ipip. The spectrum of primordial spinor density pertur- 
bations seeded during such a stage has a spectral index 
n pa 4, and is hence strongly scale dependent. In addi- 
tion, the power spectrum can be anisotropic, even though 
it is seeded within an FRW-universe. The existence of a 
"curvaton" field [15 ] during spinor-driven inflation might 
resolve these problems, ultimately resulting in the gener- 
ation of an adiabatic, nearly scale-invariant spectrum of 
density perturbations. 

Finally, for simple choices of the self-interaction V, 
there exist smooth cyclic cosmologies where the spinor 
energy density oscillates back and forth. The param- 
eters needed to accommodate a realistic cosmology do 
not appear to be unphysical. Although the simplicity of 
the models gets somewhat distorted, by a straightforward 
modification of the self-interaction V, the spinor can also 
account for dark energy and still allow for realistic cyclic 
non-singular solutions. 

Our approach has been to treat the spinor field as a 
complex valued, classical object obeying a simple rela- 
tivistic equation of motion, a non-linear generalization of 
the Dirac equation in an expanding universe. In the Ap- 
pendix we have addressed the validity and relevance of 
this assumption. Additionally, the validity of our results 
certainly depends on the stability of our homogeneous 
solutions against the growth of inhomogeneous fluctua- 
tions. We have left this question for future work. 

To conclude, we have shown that spinors can accom- 



10 



modate a large set of interesting cosmological solutions. 
Although bilinears generically decay during expansion, 
they could still be presently important if initially they 
were sufficiently displaced out of equilibrium (as for any 
interesting cosmological solution). Due to its nature, 
when dealing with a spinor field the question to ask is not 
whether certain behavior is possible, but rather, whether 
the corresponding self-interaction is natural. 
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APPENDIX A: CLASSICAL SPINORS 

A Dirac spinor is a four-component object ip that 
transforms according to (Q) and obeys Dirac's equation. 
As far as one is only concerned with solutions of equations 
of motion, the components of a spinor can be consistently 
regarded as being complex numbers, as we have done in 
this paper. However, our world is ultimately described by 
quantum-mechanical laws, and the question is: To what 
extent is a classical treatment a good approximation to 
the quantum-mechanical problem? 

In the canonical approach to quantum field theory |34j| , 
spinors are operator-valued fields that act on an appro- 
priately defined Hilbert space. The spinor operator ip 
also satisfies the Dirac equation, 

vy^d^ - mi) = 0, (Al) 

where for the purposes of illustration and simplicity we 
consider a massive fermion in flat spacetime. We work in 
the Heisenberg representation, where operators are time- 
dependent and states are time-independent. We would 
like to interpret a classical spinor as the expectation value 
of the spinor in an appropriate state |s), 

ip cl = (s\j>\s) = (i>). (A2) 

Taking the expectation value of equation (|A1|) . we find 
that ipd satisfies the equation 

i^d^d - m^ci = 0, (A3) 

which simply states that the classical spinor ip c i obeys 
the conventional Dirac equation. Therefore, we already 
recover one of the main ingredients we have used in this 
paper. 



Note that the expectation value of a spinor in a 
physical state is a complex number, not a Grassmann 
number. There exist in fact states |c) such that 
(c\ip a (x)\c) =ip a (x), with the ip a (xYs four Grassmann 
valued fields (a — 1,..,4 is the spin index). However, 
that such states are not part of the physical Fock space 
is easily seen by considering the energy density in such a 
state: 

p c = m(c\ipijj\c) — mipip . (A4) 

The energy density in a physical state must be a real 
number. However, because ip a ipb — —ipbipa, Pc=® f° r 
any n > 4. This is impossible for a non-zero real number. 

Although the expectation value of the spinor obeys the 
classical Dirac equation l|Alfl . large quantum fluctuations 
of ip around its expectation value might invalidate the 
classical approximation. In our particular case, the only 
observable that enters the classical Einstein equations is 
the energy density, which in our classical treatment is 

Pd = rmpciipd (A5) 

We want to find out whether the expectation value of 
the energy density (p) — m(ipip) is well approximated by 
(|A5I) . At this stage, one has to face a well-known prob- 
lem. The vacuum expectation value of ipip is X)fc(~ -0 = 
— oo. The conventional way of dealing with this diver- 
gence is to replace expectation values by their "renor- 
malized" counterparts, 

(...) ren = ( s |...| s )-(0|...|0). (A6) 

With this prescription, the expectation value of ipip is 
zero for the vacuum, and n for a state containing n par- 
ticles plus antiparticles per unit volume. 

Then, in order for our classical approximation to be 
valid, the following relation should hold, 

(#>ren ~ Wren Wren ^ l 
WrenW ren 

In the bosonic case, the states that satisfy inequalities 
analogous to l|A7jl have large occupation numbers. As 
the largest occupation number of fermion modes is one, 
it is commonly believed that fermionic physical states 
cannot satisfy relations such as (| AT|) . This turns out not 
to be the case. 

Let A and B be two complex numbers and let \s) be 
the state 

\s)=A\Q) + B\l). (A8) 

Here, |0) is the vacuum and |1) = Oq 1 0) is a zero- 
momentum one-particle state. The state is normalized 
if 

|A| 2 + |B| 2 = 1. (A9) 

The spinor operator can be expanded in creation and 
annihilation operators, 

ip = ^ (a-kUk + b\v k ) , (AlO) 

k 
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where Uk and Vk are normalized complex- valued spinors, 
WfeMfe = —VkVk — 1 and v^u^ = UkVk = 0. The reader 
can easily verify that for the state l)A8|) 

Wren = A* Bu , (fyren = B* Au , (^1p) r en = \B\ 2 . 

(AH) 

Therefore, condition (|A7|) implies 

\B\ 2 = 1- \A\ 2 < 1. (A12) 

Obviously, the last condition can be easily met. It means 
that a spinor can be treated classically if its quantum 
state is "close" to the vacuum. In fact, this is what one 
is doing by setting the fermions to zero in a classical 
treatment of any theory that contains spinors. But even 
if there are departures from the vacuum, we have shown 
that treating a massive Dirac spinor classically is in some 
cases a good approximation. 

Finally, let us comment on the fermion condensates 
that are often encountered in particle and condensed 
matter physics. In quantum theories with self-interacting 
fermions, it might happen that the spinor bilinear 't/jip 

9 Note, however, that one generically exp ects fermion condensates 
to couple non-minimally to gravity |3Sl . 



develops a non-zero vacuum expectation value. This is 
what occurs for instance in the BCS theory of supercon- 
ductivity j^l, where phonon- induced interactions cause 
electrons to form bound Cooper pairs. In the relativis- 
tic Nambu-Jona-Lasinio model |36| or its renormalizable 
counterpart, the Gross-Neveu model self-interacting 
chiral fermions form a scalar condensate, spontaneously 
breaking chiral-symmetry and dynamically generating a 
fermion mass. Within the effective action formalism, the 
dynamics of the condensate is completely determined by 
a classical scalar field theory. The exact classical the- 
ory that reproduces the full variety of phenomena is ex- 
tremely complicated. However, for certain states of the 
quantum system the effective action is well approximated 
by a simple, local, relativistic scalar field theory 9 . Just as 
a strongly coupled fermionic system can be effectively de- 
scribed by a classical scalar field theory, it is conceivable 
that certain strongly coupled systems might be described 
by a simple classical spinor field theory, as we consider 
in this paper. 
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